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1.0  INTRODUCTION  AND  SUMMARY 


This  work  addresses  the  possibility  of  a  high  average  power  coherent 
light  source  in  the  mid-infrared,  specifically  in  the  4.8  pm  atmospheric 
window.  The  approach  is  to  frequency-double  the  9.55  pm  CO2  laser  line, 
yielding  4.775  pm.  In  order  to  achieve  efficient  conversion,  the  intensity  of 
the  laser  is  increased  by  intracavity  resonant  modulation,  or  "mode¬ 
locking.  ’  The  intensity  for  efficient  conversion  in  typical  available  lengths 
of  doubling  material  is  at  a  minimum  20  MW  cm'^.  However,  the  most 
efficient  CO2  lasers  operate  with  a  pulse  duration  of  at  least  1  psec,  the 
characteristic  time  for  energy  transfer  from  N2  to  CO2  at  1  atmosphere. 
Without  modulation  a  fluence  of  at  least  20  J  cm'^  would  therefore  be 
needed  for  high  overall  system  efficiency.  With  modulation  at  a  1:10  or 
better  mark-to-space  ratio,  this  fluence  is  reduced  to  the  2  J  cm'^  range  that 
typifies  the  surface  damage  threshold  of  available  materials. 

One  important  unknown  was  the  effect  on  damage  threshold  of  a 
modulated  versus  an  unmodulated  pulse.  Was  the  threshold  dependent  on 
peak  intensity,  or  fluence? 

Other  questions  related  to  the  specific  material  AgGaSe2,  silver 
gallium  selenide,  chosen  for  this  work.  What  was  the  quality  of  the  second 
harmonic  beam;  what  was  the  effective  nonlinear  coefficient;  what  was  the 
temperature  dependence  of  the  phase-matching;  and  what  was  its  thermal 
conductivity  (to  remove  absorbed  heat  in  an  average  power  situation)? 

Experimentally  a  mode-locked  pulse  train  was  generated  using  a 
TEA  CO2  laser  pumped  for  a  duration  of  3  psec,  containing  1  nsec  pulses 
spaced  by  40  nsec.  The  9.55  pm  line  was  selected  either  by  injection  or  by 
the  use  of  an  intracavity  grating. 

In  a  series  of  experiments  at  increasing  fluence,  integrated 
harmonic  energy  conversions  of  up  to  78%  were  observed.  The  phase¬ 
matching  angle  was  insensitive  to  temperature  within  the  accuracy  of  our 
measurement,  the  second  harmonic  beam  was  diffraction-limited  and  the 
nonlinear  coefficient  that  best  fitted  the  data  was  dse  =  31  ±  4  pm  V  ^  in 
agreement  with  the  most  accurate  literature  measurements.  An 
interesting  observation  was  made  of  the  increased  angular  sensitivity  of 
phase-matching  at  high  harmonic  conversion. 
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The  crystal  surface  damage  threshold  was  established  for  long  and 
short  pulse  trains  of  very  different  peak  intensity.  The  damage  was 
apparently  dependent  more  strongly  on  fluence  than  intensity,  and  could  be 
represented  by 


Fq  =  Jcm“^ 

where  Fq  is  the  damage  fluence  and  t  is  the  pulse  train  duration  in 
microseconds. 

The  crystal  thermal  conductivity  was  measured  to  be  0.012  W  cm'^ 
in  directions  perpendicular,  and  mostly  parallel  to,  the  c-axis.  This  is 
a  rather  low  value,  comparable  to  that  for  fused  silica,  or  pyrex,  and  can  be 
the  cause  of  considerable  crystal  temperature  gradients  in  an  average- 
power  application. 

The  selection  of  AgGaSe2  and  a  summary  of  the  theory  of  second 
harmonic  generation  are  presented  in  Section  2.  Laser  line  selection  is 
discussed  in  Section  3.  Short  pulse  generation  by  mode-locking  is  described 
in  Section  4,  which  also  includes  a  measurement  of  the  pulse  duration  by 
autocorrelation.  The  angular  and  temperature  dependence  of  phase¬ 
matching  are  discussed  in  Section  5  and  the  detailed  comparison  of  theory 
and  experiment  are  covered  in  Section  6,  which  also  summarizes  the 
damage  measurement  and  details  the  beam  quality  measurements.  In 
Section  7  we  describe  the  thermal  conductivity  measurement.  Conclusions 
relating  to  the  use  of  AgGaSe2  for  efficient,  high  power  generation  in  the  4.8 
pm  band  are  presented  in  Section  8. 
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2.0  SELECTION  OF  AgGaSe2 


The  nonlinear  material  silver  gallium  selenide  (AgGaSe2)  was 
selected  for  these  experiments  because: 

a)  It  has  a  relatively  high  nonlinear  coefficient. 

b)  It  is  now  available  with  very  low  absorption  at  9.55  pm  (or 
10.6  pm). 

c )  It  is  available  in  perfect  single  crystals  of  length  up  to  35  mm  at 
50°  to  the  C-axis. 

d)  It  is  chemically  stable  in  room  air. 

e)  It  has  a  relatively  high  optical  damage  threshold. 

Other  nonlinear  materials  that  were  considered  are  TlsAsSes  (TAS) 
which  has  been  used  for  high  efficiency  doubling  at  10.6  pm,0-^^  and 
ZnGeP2  ‘®’^^  At  present,  there  are  still  crystal  quality  problems  in  TAS,(5) 
although  an  effective  3.3  cm  length  of  material  has  been  achieved  and  it  has 
a  very  low  absorption  at  9.55  pm.  The  second  harmonic  beam  quality  was 
not  quoted  in  References  (4)  and  (5).  ZnGeP2  has  an  absorption  in  the  10  pm 
region,  attributable  to  an  optical  phonon  overtone,  so  that  its  use  for  high 
average  power  doubling  is  questionable. 

2.1  Frequency  Doubling  in  AgGaSe2 

The  t}T)e  I  phase  matching  process^®)  is  optimal  for  frequency 
doubling  from  10.6pm  or  9.55  pm  in  AgGaSe2.  In  this  configuration  the 
fundamental  frequency  is  an  ordinary  ray  and  the  second  harmonic  is  an 
extraordinary  ray.  AgGaSe2  is  a  negative  uniaxial  crystal  with  indices 
given  by  the  following  Sellmeier  expression: 


n^  =  A  + 


B 


D 


where  n  is  the  refractive  index  and  X  is  the  wavelength  in  microns. 

The  Sellmeier  coefficients  for  AgGaSe2  are  given  by  Bhar,(9)  but  we 
Lave  found  in  the  present  work  that  they  do  not  predict  the  correct  phase¬ 
matching  angle.  The  Kildal  set  of  coefficients^!®^  has  been  modified  by 
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D.  Roberts  of  Cleveland  Crystals  to  encompass  not  only  our  present  data 
but  also  parametric  oscillator  data  from  other  laboratories.  The  Bhar  and 
Roberts  sets  are  listed  in  the  following  table. 


INDEX 

A 

B 

c 

D 

E 

no 

4.6453 

2.2057 

0.1879 

1.8377 

1600 

Bhar 

ne 

5.2912 

1.3970 

0.2845 

1.9282 

1600 

no 

3.9362 

2.9113 

0.15071 

1..7954 

1600 

Roberts 

ne 

3.3146 

3.3616 

0.14095 

1.7677 

1600 

The  type  I  phase-matching  angle  6m  for  frequency  doubling  is  given  by; 


sin^On, 


/ 

1 

f 

1 

\  ^  J 

(  \ 

1 

2 

'  \ 

1 

2 

where  subscripts  o'  and  ’e’  refer  to  ordinary  and  extraordinary  rays,  and 
subscripts  r  and  '2'  refer  to  the  fundamental  and  second  harmonic 
wavelengths.  For  doubling  of  the  CO2  P20  line  at  9.55  pm,  the  Bhar  set 
gives  0m  =  50.11°,  which  is  the  angle  to  which  our  AgGaSe2  crystal  was  cut 
(see  specification  below).  During  the  experiments  we  found  that 
0m  =  48.86  ±  0.4°  where  almost  all  of  the  error  relates  to  the  quoted  accuracy 
of  the  X-ray  crystal  orientation  performed  by  Cleveland  Crystals,  the 
supplier.  This  experimental  value  was  one  of  those  used  by  Roberts  to 
derive  the  above  modified  Sellmeier  coefficients,  which  now  give  6m  =  48.87° 
for  doubling  at  9.55  pm. 

For  later  reference,  the  following  indices  are  predicted  by  the  Roberts 
coefficients: 

ni«  =  2.596881  nie  =  2.5649089 

02"  =  2.615515  02"  =  2.582935 
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Z2  Efficdency  of  Freqiieocy  Doubling 

The  efficiency  of  doubling  is  usefully  summarized  by  Eimerl.ri2)  For 
plane  waves  the  interaction  of  the  first  and  second  harmonics  is  mediated 
by  a  coupling  constant  C  given  by: 

C  =  5.46deFF(pm^  V)  /  (A.i(|j.m)n^^^) 


where  d^PF  =  d.ig  sin  Gm  for  type  I  doubling  in  AgGaSe2  and  the  units  of  C 
are  po^  a  crystal  of  length  L,  the  conversion  efficiency  is  G3) 


n  =  tanh* 


-tanh  '(sn 


1+6^/4110 


where 

8  =  (i]dkL 

and  sn  is  a  Jacobi  elliptic  function.  The  units  of  the  first  harmonic 
intensity  1  are  GW  cm-2  for  L  in  cm.  5  is  one  half  the  phase  mismatch 
between  the  second  harmonic  waves  originating  at  the  input  and  exit  faces 
of  the  crystal. 

For  relatively  low  drive  parameters '  r|o  <  1  the  conversion  efficiency 
is; 


n  =  To 


For  larger  drive  parameters,  at  perfect  phase  matching  (5  =  0), 

Ti  =  tanh"(Tijj'-). 

The  behavior  of  r|  for  large  drives  Po  and  non- zero  6  is  discussed  by 
Eimeri,  who  gives  a  useful  contour  map  of  conversion  efficiency  in  rio,  5 
space.  The  increased  sensitivity  of  conversion  to  de-phasing  at  high  drive 
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parameters  is  discussed  below  in  Section  5.1  in  relation  to  the  present 
experiments  in  which  the  effect  was  quite  noticeable. 

2^  Doubling  with  Gaussian  Pulses 

The  mode-locked  pulses  of  the  present  experiment  have 
approximately  Gaussian  transverse  and  temporal  intensity  profiles. 
Provided  that  the  birefringent  "walk-off’  is  much  less  than  the  transverse 
size  of  the  fundamental  beam  it  is  a  good  approximation  to  integrate  the 
foregoing  plane  wave  result  across  the  range  of  intensities  at  different 
transverse  positions  in  the  beam.  For  our  doubling  process  the  estimated 
walk-off  is  0.4  mm  (in  our  35  mm  crystal)  whereas  the  1/e^  intensity  radius 
ranged  between  2  mm  and  4  mm  in  the  experiments. 

Tlie  intensity  can  be  written  as: 

where  Iq  is  the  peak  on-axis  intensity,  A  is  the  HWHM  pulse  duration  and 
w  is  the  1/e  amplitude  ( l/e^  intensity)  radius  of  the  beam. 

The  energy  in  the  pulse  is: 

oo  oe 

Ep  =  2;i  I  J  l(r,t)rdtdr 

—oo  O 

=  ioX-^^'Aw^  /  (zVInI) 

The  time  and  space  averaged  second  harmonic  conversion  efficiency  is,  for 
5  =  0: 


n  =  ^  7  T  2jirl(r.t)tanh~(riJ,''2)dtdr 

tp  -oo  O 

This  function  is  shown  vs  Figure  1  (labelled  "Gaussian").  Also 

shown  is  the  value  of  r|  for  the  flat-topped  plane  wave  case  discussed  above 
in  Section  2,2.  This  figure  shows  that,  for  example,  to  obtain  50%  energy 
conversion  with  a  Gaussian  pulse  requires  about  twice  the  drive  parameter 


-6- 


PEAK  SHG  DRIVE  PARAMETER  CL 

P24J5 


Fig.  1  Second  harmonic  conversion  efficiency  as  a  function  of  the  drive 
parameter  Vrio  =  CL  VI,,  where  Iq  is  the  peak  intensity  in  a  spatial  and 
temporal  Gaussian  profile.  Dashed  line  shows  constant  plane  wave  case. 
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rioL/2  and  about  four  times  the  intensity  needed  for  the  same  conversion 
from  a  hypothetical  flat-topped  plane  wave  pulse. 

Later,  in  comparing  the  measured  conversion  efficiency  of  pulse 
trains  with  theory,  a  range  of  different  peak  pulse  intensities  (Iq)  has  to  be 
summed. 

2.4  Anticipated  Nonlinear  Coefficient  of  AgGaSe2 

For  AgGaSe2  there  is  a  range  of  dss  values  in  the  literature.  Values 
obtained  by  direct  methods,  and  also  one  comparison  measurement  are 
listed  in  the  following  table  which  includes  our  present  result: 


dselpmV-i) 

METHOD 

REFERENCE 

33  ±4 

Doubling  lO.fiiim  in 
comparison  with  known 
d36  for  GaAs 

(2) 

32.4  ±5 

Doubling  10.6|i,  direct 

measurement 

(10) 

43  ±(?) 

Doubling  10.6M,m,  direct 

measurement 

(14) 

39  ±  (?) 

Parametric  gain  at 
3.39pm,direct 

measurement 

(15) 

31  ±4 

Doubling  9.55|j.m, 
direct  measurement 

present  work 

The  details  of  our  present  measurements  are  discussed  in 
Section  6.0,where  the  data  is  fitted. 


2.5  Ciystal  Specification 

The  following  specification  was  agreed  with  Cleveland  Crystals,  in 
which  0  is  the  polar  phase-matching  angle  (Bhar  Sellmeier  used)  and  ip  is 
the  azimuthal  phase-matching  angle. 
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DESCRIPTION 

1  AgGaSo2  -  10  X  10  X  35  mm  - 1  -  0  =  45°,  0  =  50.11° 

SHG  -  9.550/4.775  ^im 

Single  Crystal 

Polish  (80%  of  aperture): 

Flatness  -  best  efforts 

Parallelism  -  ±  15  arc-seconds 

Surface  Roughness  ~  25  A  rms  -  best  efforts 

AR  "W"  coated  both  surfaces  for  9.550  and  4.775  pm 

No  major  inclusions  and  best  effort  on  minimum  of 

minor  ’nclusions  (major  inclusion  size  is  greater  than 

50  pm) 

1  Calonmetnc  Absorption  Measurement  at  10.6  pni 

1  Profilometry  (Surface  Roughness) 

1  IR  Photography  (NIR  Broadband,  B&W)  showing 

propagation  direction  and  also  a  photo  of  the  slab 
before  cutting  to  five  the  location  of  the  crystal  in  the 
boule  and  venfy  its  single  crystal  nature 

Shipping  Schedule:  6  Weeks  AKO 

The  crystal  was  delivered  on  schedule  accompanied  by  the  following 

data; 

1.  Infrared  photograph  at  1  pm  showing  absence  of  inclusions 

(Figure  2  (a)). 

2.  Sketch  showing  location  of  crystal  cut  in  boule  (Figure  2  (b)). 

3.  Surface  Profiles  at  a  number  of  locations  in  the  central  80%  of  the 

crystal  aperture  (Figure  3),  showing  approx.  20  A  r.m.s.  finish. 

The  absorption  at  10.6  pm  was  measured  to  be  0.012  cm'^  by 

Cleveland  Crystals.  In  order  to  evaluate  the  "W"  anti-reflection  coatings, 

we  ran  the  crystal  on  a  Perkin  Elmer  infrared  spectrophotometer.  At  9.55 

pm  and  at  4.775  pm  the  transmission  was  90%  ±  2%.  An  additional 
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SURFACE  PROFILE 


SURFACE  PROFILE 


measurement  with  the  polarized  C.W.  laser  at  9.55  pm  gave  a  transmission 
of  92.5%  ±  1%.  Considering  the  absorption  loss  of  4%  at  10.6  pm,  which  is 
probably  the  same  at  9.55  pm  because  residual  "absorption "  in  AgGaSe2  is 
spectrally  featureless  C6)  (possibly  due  to  scattering  by  fine  inclusions),  the 
coatings  were  contributing  a  loss  of  about  2%  at  each  surface  at  9.55  pm. 
Because  scattering  increases  at  reduced  wavelength,  it  is  likely  that  the 
absorption"  was  at  least  4%  at  4.775  pm  and  that  the  coatings  each 
reflected  less  than  3%  at  this  wavelength.  The  coatings  themselves 
consisted  of  a  multilayer  combination  of  lead  and  thorium  fluorides.  In 
Section  6.2,  the  optical  damage  threshold  of  the  coatings  is  given  for  various 
incident  pulse  train  profiles. 
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3.0  SELECTION  OF  9.55  LINE 


The  selection  of  the  9.55  Mjn  P20  CO2  laser  line  was  performed  by 
either  line  injection  of  a  CW  laser,  or  by  grating  tuning.  Both  techniques 
yielded  comparable  mode-locked  pulse  durations  but  the  grating  technique 
proved  simpler  and  more  reliable.  Injection  had  first  been  considered  as  a 
technique  for  eliminating  the  gain-switched  "spike"  in  order  to  obtain 
multiple-microsecond  pulse  trains,  for  harmonic  generation,  of  relatively 
uniform  pulse  intensity.  Although  this  was  achieved  with  very  careful 
tuning, '  and  a  feedback  monitor  to  establish  resonance,  the  practical 
usefulness  of  the  approach  was  not  adequate.  By  changing  the  gas 
mixture,  carefully  setting  the  cavity  gain  and  amplifying  the  "tail"  only, 
satisfactory  long  pulse  waveforms  were  obtained  with  much  better 
reliability  from  the  grating-tuned  laser. 


3.1  Lqjection-Locking 

The  arrangement  for  injection-locking  is  shown  in  Figure  4.  A  CW 
CO2  laser  gave  an  output  of  2  W  in  a  TEMoo  mode  at  9.55  pm.  The  CW  beam 
was  allowed  to  expand  for  a  distance  of  427  cm  before  passage  through  a  500 
cm  collimating  lens.  The  beam  emerged  parallel,  with  a  calculated  l/e^ 
spot  parameter  of  0.29  cm,  slightly  smadler  than  the  0.46  cm  oscillator  spot 
parameter  at  the  output  mirror  of  the  oscillator,  where  injection  occurred. 

A  957c  transmitting  beam  splitter  was  used  to  direct  0.2  W  of  CW  power  into 
the  oscillator  (or  at  least  in  the  direction  of  the  oscillator).  The  same  beam 
splitter  transmitted  95%  of  the  oscillator  output  for  use  in  the  experiment, 
and  although  5%  of  the  output  was  focussed  back  to  the  CW  laser,  the 
fiuence  there  was  never  more  than  0.2  J  cm*^,  so  the  CW  laser  optics  did  not 
suffer  damage. 

The  P20  9.55  pm  line  has  less  gain  than  about  10  of  the  10.4  pm  band 
lines  and  several  other  9.4  pm  band  lines.^^^^  Without  line  selection  the 
CO2  laser  emits  10.6  pm  -  a  line  with  a  gain  of  1.15  times  that  at  9.55  pm. 
Gain  on  the  10.4  pm  band  could  be  reduced  by  the  use  of  a  thin  CaF2 
window  at  Brewster  angle  within  the  oscillator  cavity.  The  0.92  mm  thick 
window  at  incidence  angle  54.3°  transmitted  an  estimated  64%  at  10.6  pm, 
and  91%  at  9.55  pm,  completely  suppressing  all  10.4  pm  band  emission. 
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OUTPUT  TO 
AMPUriER 

P2W0  OR  CRYSTAL 


Fig.  4  Optical  arrangement  for  the  injection-locking  of  the  CO2  oscillator 
at  9.55  pm. 


-14- 


With  the  CaF2  window  but  no  other  line  selection,  emission  was  mostly  on 
the  P18  and  P16  lines  of  the  9.4  (im  band,  with  only  occasional  emission  on 
the  desired  P20  line. 

Line  selection  via  injection  was  found  to  depend  on  careful  angular 
alignment  and  on  fine  tuning  of  the  oscillator  cavity  length  to  within  a 
fraction  of  the  CO2  laser  wavelength.  The  latter  adjustment  was  performed 
manually  \ia  a  fine  linear  drive  which  carried  the  laser  output  mirror.  A 
problem  was  that  this  adjustment  would  only  persist  for  at  best  a  few  tens  of 
seconds,  before  either  vibrations  or  thermal  changes  altered  the  cavity 
length.  In  order  to  achieve  predictable  results  a  monitor  of  the  intracavity 
injected  intensity  was  set  up,  consisting  of  a  chopper  and  pjn’oelectric 
detector  positioned  to  intercept  the  reflection  off  one  face  of  the  modulator 
(slightly  tilted  to  the  cavity  axis).  When  this  monitor  signal  had  been 
optimized  by  fine  tuning  of  the  cavity  length,  the  oscillator  pulse  was  always 
line-locked.  Varying  degrees  of  optimization  were  employed  to  generate  the 
range  of  pulse  waveforms  illustrated  in  Figure  5.  The  ratio  of  "spike"  to 
"tail"  intensities  could  be  reduced  to  less  than  3:1  in  the  best  case.  In  this 
case  the  delay  to  oscillator  output  was  reduced  to  2.0  psec,  whereas  with  no 
injection  this  delay  lengthened  to  3.0  psec  (and  the  9.55  pm  energy  reduced 
almost  to  zero). 

3J2  Line  Selection  by  Grating 

A  standard  high  power  CO2  laser  grating  was  obtained  from  Milton 
Roy  of  Rochester,  NY.  It  was  gold-coated,  had  100  grooves  per  mm  and  was 
blazed  for  9.5  pm,  at  which  wavelength  its  reflectivity  was  nominally 
greater  than  95%.  The  incidence  angle  for  9.5  pm  is  25.5°. 

Neighboring  CO2  laser  lines,  P18  and  P22,  are  reflected  with  an 
angular  deviation  of  0.04°  which  translates  after  a  single  600  cm  roimd  trip 
to  a  lateral  displacement  of  0.4  cm,  which  equals  nearly  half  of  the  mode- 
selection  aperture  diameter  and  therefore  gives  heavy  losses.  Over  the 
3  GHz  bandwidth  of  the  P20  line  itself  the  displacement  after  one  round  trip 
is  only  0.02  cm,  so  that  about  10  round  trips  would  be  required  to  introduce 
only  a  few  percent  of  loss.  The  grating  therefore  does  not  significantly 
narrow  the  laser  spectrum  within  a  given  line. 

As  expected,  the  gain-switched  "spike"  dominates  the  laser  output 
( for  the  3  psec  PFN  that  we  are  using)  as  shown  in  Figure  6  (a).  By 
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Fig.  5  9.55  (im  output  intensity  profiles  of  injected  laser  with  various 

degrees  of  cavity  length  matching;  -  (a)  best  length  match;  (b)  (d) 

decreasing  quality  of  match;  (e)  no  injection. 
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Fig.  6  Different  mode-locked  pulse  train  envelopes  achieved  with 
combinations  of  a  grating-tuned  oscillator  and  amplifier:  -  (a)  maximum 
oscillator  aperture  consistent  with  TEMqo  mode;  (b)  reduced  oscillator 
aperture;  (c)  oscillator  output  with  enhanced  N2  in  mix;  (d)  output  (c) 
amplified  with  a  2.5  |isec  delay. 


reducing  the  mode-selection  aperture,  the  cavity  losses  could  be  increased 
to  the  extent  that  lasing  ceased  immediately  after  the  spike  (Figure  6  (b)) 
leaving  a  useful  short  pulse  train  for  harmonic  generation  at  high  intensity 
and  low  fluence.  By  increasing  the  nitrogen  content  of  the  mix  (from  our 
standard  He:C02:N2  =  8:2:1  to  7:2:2)  the  "tail"  could  be  enhanced  relative  to 
the  spike  (Figure  6  (c));  and  following  amplification  with  a  2.5  psec  delay  on 
the  amplifier  trigger  pulse,  the  train  could  be  dominantly  "tail"  as  shown 
in  Figure  6  (d).  The  latter  pulse  train  was  useful  for  the  demonstration  of 
harmonic  generation  at  moderate  intensity  and  high  fluence. 
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4.0  MODE-LOCKING 


The  unmodulated  pulsed  CO2  laser  emits  a  noisy  pulse  train,  in 
which  there  is  periodicity  at  the  cavity  round  trip  time  period.  Examples  for 
the  present  oscillator  are  shown  in  Figures  7(a)  and  7  (b).  The  duration  of 
individual  "spikes  "  within  the  train  is  approximately  2  nsec.  This  number 
is  the  result  of  gain  narrowing  of  the  laser  spectrum  during  the  pulse 
buildup  from  noise.  The  spontaneous  emission  bandwidth  of  a  1 
atmosphere  CO2  laser  is  about  3  GHz,  so  that  an  initial  noise  pulse  duration 
of  0.1  nsec  must  have  existed. 

We  employed  amplitude  modulation  using  diffraction  off  an  acoustic 
wave  in  a  germanium  crystal.  The  availability  of  an  oscillator,  amplifier 
and  several  Ge  crystals  with  transducers  at  12.5  MHz  dictated  our  choice  of 
cavity  round  trip  time  (40  nsec)  and  modulation  frequency.  The  acoustic 
wave  amplitude  falls  to  zero  twice  in  each  cycle,  so  that  an  effective 
25  MHz  (=  40  nsec)  frequency  is  achieved. 


4.1  Modulator  Test 

The  modulators  were  characterized  in  a  separate  benchtop  test  in 
which  a  3W  CW  9.55M.m  beam  was  passed  through  the  modulator  and  then 
a  small  iris  at  3  m  beyond  the  modulator.  The  transmitted  intensity 
decreased  as  a  function  of  RF  power  to  the  modulator  as  shown  in  Figure  8. 

The  diffraction  efficiency"’  is  given  by; 


where 

II 

lo 

Pa  = 
M2(Ge  = 
L 
X 


_2  1 2 

— M2PA  ^ 
2  ^ 


1/2 


diffracted  power 
incident  power 

acoustic  power  density  (ergs  cm'2  sec*^) 

8.1  X  10' 16  sec^  g-i 

length  (cm)  through  Ge  crystal  (optical  direction) 
wavelength  of  incoming  light  (cm) 
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Fig.  7  Spontaneous  modulation  of  the  grating-tuned  CO2  laser  output 
intensity.  40  nsec  periodicity  is  seen  in  each  example. 
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Fig.  8  Transmission  of  the  undeflected  component  of  a  9.55  pm  CW  laser 
as  a  function  of  modulator  RF  drive  power. 
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The  above  formula  is  correct  in  the  "Bragg  regime  ’  for  which: 


where  A  =  acoustic  wavelength.  For  Ge  the  speed  of  sound  is 
5.5  X  10^  cm  sec"^  in  the  111  direction  which  has  the  highest  acousto-optical 
figure  of  merit  M2  for  light  polarized  in  the  111  direction.^  At  12.5  MHz 
the  acoustic  wavelength  is  therefore  A  =  0.044  cm.  The  available  crystals 
had  L  =  1  cm,  hence; 

A-^ 

and  the  Bragg  regime  applied. 

The  benchtop  modulator  test  measured  the  time-averaged  diffraction 
efficiency  due  to  an  acoustic  standing  wave.  In  such  a  case  the  diffraction 
efficiency  is  periodic  in  time  and  an  average  has  to  be  taken.  The  acoustic 
power  density  Pa  is  proportional  to  the  square  of  the  acoustic  amplitude, 
and  the  standing  wave  amplitude  goes  as  sinlwgt)  where  Wg  is  the  acoustic 
angular  frequency.  Therefore,  the  functional  form  of  the  diffraction 
efficiency  with  time  is: 


|j-  =  sin“[kASin(Wst)] 


Averaging  over  one  half  acoustic  period,  and  changing  variables  to 
y  =  kA  sinwst; 


-Ws  ’T"  l|  ^  2  V  sin-y  ^ 

— S-  I  -^dt  =  -  J  - - ^dy 

^  o  'o  ^  0  yl  —  y~ 


From  the  modulator  test,  for  example,  a  time-averaged  diffraction 
efficiency  of  50%  is  observed  at  a  modulator  power  of  26  Watts,  and  this 
corresponds  to  kA  =  0.55. 
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4^  Theoretical  Elstiinate  of  Mode-Locked  Pulse  Ehiratioii 

From  the  model  of  Kuizenga  et  the  mode-locked  laser  emits 

pulses  of  (FWHM)  duration: 

f7n2  1  1 

^  2  Jt 

where 

M  =  number  of  round  trips  during  "buildup " 

0  m  =  depth  of  modulation 

fm  =  modulator  frequency  (here  fm  =  12.5  MHz) 

This  result  applies  to  the  "transient "  case,  before  a  bandwidth- 
limited  pulse  is  formed,  and  it  can  be  shown  that  the  present  work  is  in  this 
regime. 

The  single  pass  modulator  transmission  assumed  by  Kuizenga  is: 


f(t)  =  cos“(0mSinWmt). 

For  an  acoustic  standing  wave  we  observed  (above)  a  deflected  intensity: 

I,  =loSin"(kASinWst), 


or  a  transmitted  intensity  of; 

I,  =loCOS^(kASinwst). 

implying  the  identity  0m  =  kA  for  the  modulation  depth.  kA  is  already 
known  from  the  above  modulator  test  (Section  4.1). 

The  duration  of  individual  mode-locked  pulses  was  measured  using 
fast  photodetectors  at  different  points  in  the  pulse  train  and  at  different 
deptiis  of  modulation.  In  the  case  of  short  pulse  trains,  consisting  of  only 
about  10  pulses,  in  the  gain-switched  spike,  the  duration  was  also 
measured  in  a  second  harmonic  autocorrelation  experiment.  Overall,  a 
rough  agreement  was  observed  between  measured  and  theoretically- 
predicted  pulse  durations.  The  data  is  summarized  and  compared  with 
theory  in  Table  1.  The  measured  pulse  duration  was  always  less  than  the 
predicted,  by  a  factor  of  0.6  -  0.7,  for  reasons  not  understood  at  present. 
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Table  1 

Comparison  of  Predicted  and  Measured  Pulse  Durations 
for  Amplitude-Modulated  CO2  Laser 


Depth  of 
Modulation 

0m 

Number  of  Round 
Trips  from  Start 
of  Gain  (M) 

Predicted  Pulse 
Duration  tp  (nsec) 

Measured  Pulse 

Duration  (nsoc) 

0.55 

50 

3.2 

2.4 

0.55 

100 

2.3 

1.7 

0.60 

50 

2.9 

1.9 

0.60 

100 

2.1 

1.4 

1.0 

60 

1.6 

0.9-1.0 

4^  Measurements  of  Mode-Locked  Pulse  Duration 
A)  Using  Photodetector 
This  set-up  involved: 

1.  A  fast  magnetically  biased  HgCdTe  detector  (Boston  Electronics 
Model  P005-4)  with  a  bandwidth  of  500-800  MHz. 

2.  60  cm  length  of  SMA  cable  rated  at  18  GHz. 

3.  Tektronix  7A29  pre-amplifier  rated  internally  at  1000  MHz  (0.35  nsec 
risetime)  and  jielding  750  MHz  (0.47  nsec  risetime)  in  the  7104 
mainframe. 

4.  Polaroid  photograph  at  sweep  speed  of  2  ns/div. 

The  anticipated  risetime  of  this  arrangement  is  approximately  0.5  nsec.  In 
practice,  pulses  with  a  full  width  at  half  maximum  (as  photographed)  of 
0.95  nsec  were  the  shortest  observed,  under  optimum  conditions  of  mode¬ 
locking  (Figure  9).  This  allows  us  to  place  an  upper  bound  on  the  system 
risetime  of  0.475  nsec.  Most  of  the  frequency  doubling  measurements  were 
made  with  a  measured  (FWHM)  pulse  width  of  1.3  nsec,  which  is 
deconvoluted,  using  0.475  nsec  in  quadrature,  to  1.03  nsec.  Because 
1.3  nsec  is  not  clearly  larger  than  the  system  response  time,  there  is  quite  a 
lot  of  uncertainty  attached  to  this  deconvolution.  A  second,  more  accurate. 
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Fig.  9  Oscilloscope  traces  of  individual  short  pulses. 
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method  of  measurement  involves  autocorrelation  using  second  harmonic 
generation,  discussed  in  the  following  section. 

b )  rising Autocorrrelation 

In  this  technique  the  mode-locked  CO2  laser  beam  is  split  into  two 
equal  beamlets,  one  of  which  is  delayed  by  varying  amoimts  relative  to  the 
other  before  they  are  both  superimposed  within  a  nonlinear  crystal.  We 
chose  a  set-up  in  which  both  beams  had  the  same  polarization  in  order  to 
utilize  the  35  mm  AgGaSe2  crystal  in  the  Type  I  phase  matching 
configuration.  This  has  the  minor  disadvantage  that  when  one  pulse  is 
delayed  by  more  than  a  whole  pulse  length  from  the  other,  there  is  still  a 
background  level  of  harmonic.  Also,  the  two  beams  had  to  be  tilted  at  a 
small  angle  to  each  other  in  order  to  avoid  Michelson  fringes  which  would 
otherwise  give  a  fluctuating  harmonic  signal  depending  on  the  number  of 
wavelengths  of  delay.  We  used  an  intersection  angle  of  0.005  rad,  which 
>Telded  3  fringes  across  the  combined  beam,  in  the  horizontal  direction, 
perpendicular  to  the  polarization.  The  transverse  distance  between  fringes 
was  2  mm,  which  exceeded  the  calculated  "walk-off’  distance  (for  SHG)  of 
0.35  mm,  so  that  the  whole  length  of  the  crystal  was  efficiency  used.  The 
total  angle  was  less  than  half  the  available  phase  matching  angle. 

Harmonic  conversion  efficiency  was  kept  below  15%  so  as  to  be  in  the  second 
power  law  regime.  The  expected  autocorrelation  waveform  for  Gaussian 
pulses  is  now  derived. 

Where  A  is  the  HWHM  of  the  Guassian  pulses,  the  first  and  second 
pulses  have  the  intensity  waveforms: 

M  -  6  T 

A" 

A^ 

where  I2  is  delayed  by  time  t  with  respect  to  Ii.  The  wave  amplitudes  are 

E,  =  ^  and  E2  =  ^ 

(in  magnitude). 
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For  beams  superposed  at  a  small  angle  the  resxiltant  amplitude  is 
(El  +  E2  sin  <1) )  for  -  oo  <  t  <  1/2  and  {E2  +  Ei  sin  0)  for  x/2  <  t  <  «»  where  0  is  a 
linear  function  of  the  transverse  coordinate  perpendicvilar  to  the  lines  of  the 
fringes.  It  may  be  seen  from  Figure  10  (a)  that  there  is  time  symmetry 
about  t/2,  so  the  total  second  harmonic  energy  is: 

=  {E,+E2Sin0rd0dt 

0  -00 


where  the  average  is  taken  over  one  whole  fringe  in  the  transverse 
direction.  Expanding  the  bracket  and  integrating  in  0  gives: 


Esh  =  2  1  e  > 


dt 


dt 


t/2 

+0.75  I  e- 


-2/n2(t-x)^ 


dt 


This  expression  can  be  evaluated  numerically  as  a  function  of  x/A 
(Figure  10  (b))  and  compared  with  the  measured  autocorrelation  profile  to 
fit  A  (Figure  11).  The  limiting  case  of  widely  separated  pulses  (x  »  A)  gives 


^SH  = 


which  of  course  equals  twice  the  second  harmonic  energy  from  an 
individual  pulse  '  2  “i  Ej'dtl. 
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normalised 


delay 


Fig.  10  (a)  Derivation  of  autocorrelation  integral  using  symmetry  about 
time  t/2.  (b)  Computed  harmonic  signal  as  a  function  of  normalized  delay 
x/A  where  A  is  the  HWHM  of  the  Gaussian  pulse. 


harmonic  by  a  factor  of  2.1875. 

This  theory  assumes  that  all  of  the  emitted  laser  spectrum  is  mode- 
locked,  so  that  there  is  temporal  coherence  throughout  the  laser  pulse,  for 
the  leading  and  trailing  parts  of  the  pulse  to  interfere  in  a  fringe  pattern, 
regardless  of  the  delay.  The  observed  ftinge  pattern  at  most  delays  has  a 
high  contrast  ratio  (on  heat  sensitive  paper,  or  liquid  crystal  polymer)  so  the 
condition  appears  to  be  satisfied.  Additionally,  the  frequency  spectrum  of 
the  mode-locked  train  is  estimated  to  be  a  substantial  fraction  of  the 
available  gain  bandwidth  of  the  1  atmosphere  CO2  laser  medium,  so  that 
unlocked  modes  are  not  likely  to  be  present. 

The  experimental  autocorrelation  profiles  fit  the  theory  reasonably 
well  (Figure  11).  Three  determinations  of  the  peak  to  background  ratio  gave 
an  average  of  2.2  ±  0.1,  in  agreement  with  theory. 

The  measured  pulse  durations  varied  somewhat  after  different 
tunings  of  the  laser,  and  averaged  0.90  ±  0.05  nsec,  for  the  short  pulse 
trains  used  in  a  detailed  comparison  of  second  harmonic  with  theory 
(Section  6.0). 


4.4  Mode-Locking  with  Line  Selection 

Two  methods  were  used  to  select  the  P-20  9.55  pm  line  chosen  for  this 
work,  injection  by  CW  laser  and  intracavity  grating  tuning.  A  description 
of  each  set-up  is  given  in  Section  3.  In  this  section  the  suitability  of  each  of 
these  approaches  for  use  with  an  intracavity  modulator  is  discussed.  The 
pulse  durations  measured  by  each  approach  were  in  fact  indistinguishable, 
for  comparable  depth  of  modulation.  The  difference  between  them 
theoretically  is  that  in  the  line  injection  case  the  initial  intracavity  flux  has 
a  much  smaller  bandwidth,  expected  to  be  of  the  order  of  1  MHz,  whereas  in 
the  grating  case  the  initial  (noise)  bandwidth  is  on  the  order  of  1  GHz.  In 
each  case  the  modulator  is  powered  200  psec  before  the  CO2  gain  medium  is 
pulsed,  so  that  in  theory  the  intracavity  line-injected  radiation  has  some 
modulation  prior  to  the  appearance  of  gain.  This  effect  is  not  very  large, 
though,  because  the  lifetime  of  cavity  photons  without  gain  is  only  a  few 
round  trips,  compared  to  the  eventual  persistence  of  50-100  round  trips 
during  the  gain-switched  buildup  phase. 
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The  2  nsec  structure  present  in  an  immodulated  grating-tuned  laser 
(Figure  7)  is  the  result  of  gain-narrowing  during  buildup.  Most  of  this 
narrowing  is  due  to  the  gain  profile  of  the  laser  medium,  but  some  of  it  is 
due  to  grating  dispersion.  Without  this  gain-narrowing  the  individual 
pulses  wotdd  emerge  from  a  modulated  laser  with  a  1-2  nsec  envelope  (as  in 
the  line  injection  case)  but  there  would  be  finer  structure  (0.1  nsec)  within 
the  envelope.  For  the  CO2  laser,  the  relatively  narrow  gain  bandwidth 
ensures  that  with  the  present  frequency  and  depth  of  modulation  the 
resultant  pulses  are  unstructured  and  of  similar  duration  in  both  the  line 
injection  and  grating  timed  cases.  This  need  not  necessarily  have  been  the 
case.  For  example,  the  line-injected  and  mode-locked  Nd:YAG  laser  would 
not  have  structured  pulses,  but  structure  would  usually  be  present  in  the 
grating-tuned  version.  In  the  transient  mode-locking  regime  the  duration 
of  the  pulse  envelope  would  be  much  the  same  in  each  case. 

One  aspect  relating  to  the  use  of  a  grating  within  a  modulated  cavity 
is  that  apparently  the  round  trip  time  is  different  at  different  lateral 
positions  on  the  grating  (t)T)ically  used  at  about  25°  to  normal).  For  a  spot 
size  of  0.6  cm,  the  difference  in  cavity  roimd  trip  length  is  about  1  cm.  First 
of  all,  the  mode-locked  pulse  duration  was  experimentally  insensitive  to 
±  0.5  cm  variation  in  cavity  length  (out  of  a  round  trip  length  of  600  cm),  but 
more  importantly  the  modes  of  a  cavity  with  a  grating  are  more  like 
supermodes "  in  a  distributed  feedback  resonator.  The  round  trip  length 
becomes  "fuzzy  ”  and  the  only  effect  to  be  considered  is  the  small  difference 
in  cavity  loss  across  the  CO2  laser  line  due  to  misalignment,  because  the 
grating  only  reflects  "normally"  at  the  central  frequency. 

4.5  Modulator  Damage  Threshold 

The  1  cm  thick  Ge  modulator  crystals  occasionally  exhibited  surface 
optical  damage  (probably  to  coatings)  and  at  other  times  laser  energy  was 
down  but  could  be  restored  by  moving  the  modulator  transversely  to  an 
unused  position  indicating  possible  volume  damage,  even  although  surface 
damage  was  not  always  visible. 

Through  trial  and  error  it  was  found  that  hundreds  to  thousands  of 
pulses  could  be  obtained  without  apparent  damage  for  intracavity  fluences 
on  the  modulator  of  up  to  0.45  J  cm‘2.  This  result  refers  to  mode-locked 
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pulse  trains  of  10-100  pulse  duration,  at  peak  intensities  of  50-20  MW  cm'^, 
respectively. 

The  interpretation  of  these  numbers  is  complicated  by  the  fact  that 
the  modulator  is  very  close  to  the  grating  at  one  end  of  the  cavity  and  the 
1  nsec  pulses  simultaneously  pass  through  the  modulator  in  both 
directions,  each  way  having  a  fluence  of  0.5  J  cm'2.  The  bench  test  of  a 
modulator  should  be  a  single  pass  fluence  of  1  J  cm'2.  Further  uncertainty 
comes  from  the  likely  existence  of  long  duration  standing  waves  in  which 
the  electric  field  maximxim  dwells  locally  in  parts  of  the  modulator. 

Without  anti- reflection  coatings  the  Ge  modulator  would  have 
(incoherent)  surface  losses  of  36%  per  surface.  Use  of  Brewster  angle 
modulators  is  not  practical  because  0b  =  76°  and  too  large  an  optical 
footprint  would  result. 

Damage  to  other  elements  in  the  optical  cavity  was  only  observed 
above  the  modulator  threshold.  The  ZnSe  Brewster  angle  windows  of  the 
oscillator  withstood  fluences  of  1.5  J  cm*2  for  hundreds  to  thousands  of 
pulses.  The  grating  showed  what  appeared  to  be  a  small  amount  of 
additional  finely-distributed  scattering  in  the  area  of  irradiation  at  fluences 
above  about  1  J  cm'2,  but  this  did  not  appear  to  diminish  its  efficiency,  and 
did  not  progress  into  more  advanced  damage. 
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5.0  ANGULAR  AND  TEMPERATURE  DEPENDENCE 
OF  PHASE-MATCHING 


In  a  high  average  power  application  residual  absorption  can  raise 
the  temperature  of  the  crystal,  leading  to  refractive  index  changes  which 
usually  cause  a  loss  of  phase  matching.  There  had  been  two  reports  (20,15) 
of  widely  differing  temperature  dependence  of  the  index  difference  in 
AgGaSe2,  so  to  resolve  this  conflict  the  low-intensity  phase  matching  profile 
was  established  as  a  function  of  temperature  for  the  doubling  of  mode- 
locked  pulse  trains  at  9.55pm.  This  work  was  done  with  "hybrid”  pulse 
trains  consisting  of  a  "spike"  and  "tail"  of  roughly  equal  energy,  as 
illustrated  in  Figure  6.  The  harmonic  conversion  efficiency  was  about  5% 
for  each  of  the  experimental  runs  shown  in  Figure  12.  The  doubling  was 
performed  with  a  TEMoo  beam  of  1/e^  (intensity)  radius  =  0.50  cm,  and 
almost  plane  wave  front.  The  angular  dependence  of  the  doubling  was  due 
to  the  induced  phase  mismatch  and  not  to  angular  components  in  the  first 
harmonic  beam. 

In  Figure  12  the  continuous  curves  represent  a  least  squares  fit  of  the 
(sin  6/5)2  function  to  the  data  points,  with  the  appropriate  scaling  and  offset 
parameters.  Here 


5  = 


Iakl 

2 


where  L  =  crystal  length  and 

Ak  =  k"«-2k'^=^{n®-n°) 

where  k  =  9.55  pm. 

In  four  fitting  procedures  the  mean  value  of 

^(nf-nO)e.e„ 

was  established  as  -0.0325  ±  0.0025  per  sterad,  with  the  sign  appropriate  for 
0  increasing  with  increasing  angle  to  the  optic  axis,  6  being  measured 
internal  to  the  crystal. 
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The  same  quantity  may  be  estimated  theoretically,  starting  from  the 
relation 


nf(e) 


'"cose'' 


■2  j 


f  . 

Sine 


V  “2  ) 


and  differentiating  with  respect  to  e 


Since 


0n° 


-2  _  -2  cose  sine  2  sine  cose 

^  Hf 


is  zero, 


— k-n°) 

ae^  -  '  / 


Theor 


-0.0322 


which  is  equal  to  the  measured  value  above  within  the  experimental  error. 
For  the  record,  the  index  for  both  "e”  and  "o"  rays  at  the  phase  matching 
angle  is  2.59693. 

As  may  be  seen  in  Figure  12,  there  is  not  any  detectable  change  in  the 
phase  matching  angle  with  temperature  between  30°C  and  70°C.  In  part 
(a)  of  that  figure  the  matching  is  0.025°  (external)  further  from  the  C-axis  at 
the  higher  temperature,  with  a  fitting  error  of  ±  0.02°.  In  part  (b)  the 
matching  is  0.03°  ±  0.02°  closer  to  the  C-axis  at  the  higher  temperature. 

The  average  of  these  measurements  is 

^(nl-n'>)e.e„=0±2xl0-^K-> 

The  result  is  much  lower  than  two  previous  determinations,  which 
however  relate  to  different  wavelengths  or  temperature  ranges.  In 
reference  (1^)^  a  value  of  1.1  x  10  ®  K‘l  was  determined  for  doubling  of 
10.6  pm  between  lOOK  and  298K,  whereas  in  reference  (20)  a  single 
wavelength  of  3.39  pm  was  used,  the  temperature  range  was  lOOK  to  300K 
and  a  value  of  30  x  10‘®  K'l  was  obtained. 
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To  conclude  this  section,  the  angular  dependence  of  phase  matching 
was  exactly  equal  to  that  theoretically  expected  for  a  3.5  cm  crystal, 
indicating  the  absence  of  index  flaws  throughout  its  length.  Furthermore, 
the  temperature  coefficient  of  the  index  difference  was  zero  to  within  the 
experimental  error,  suggesting  that  we  have  accidentally  hit  upon  a  null  at 
which  the  temperature  coefficients  of  ordinary  and  extraordinary  indices 
are  exactly  equal.  Individually  these  coefficients  are  roughly  in  the  range 
40-70  X  10"®  K'l  at  room  temperature.^ 

5.1  Intensity  Dependence  of  Phase-Matching 

It  is  known  that  as  the  conversion  efficiency  into  second  harmonic  is 
increased  the  process  becomes  more  sensitive  to  small  degrees  of  phase 
mismatch.^21,12)  This  effect  could  readily  be  observed  in  the  present 
experiments.  Figure  13  contrasts  the  angular  dependence  of  second 
harmonic  efficiency  at  low  conversion  (about  5%)  with  that  at  high 
conversion  (68%  energy  conversion  at  zero  phase  mismatch).  The  width  of 
the  central  peak  (at  half  height)  has  been  reduced  by  a  factor  of  0.67  in  this 
run,  which  refers  to  the  longer  pulse  train  shown  in  Figure  6,  at  an 
incident  energy  of  140  mJ.  The  secondary  maxima  are  higher  and  closer  to 
the  central  peak,  in  line  with  expectation.  Another  experimental  run 
showing  a  width  reduction  by  a  factor  of  0.57  is  shown  in  Figure  14.  This 
run  had  a  smaller  energy  conversion  efficiency  at  zero  phase  mismatch 
(56%)  and  employed  the  short  pulse  trains  shown  in  Figure  6  at  an  incident 
energy  of  70  mJ.  The  small  differences  between  the  two  data  sets  could  be 
due  to  deviations  from  Gaussian  transverse  profiles  (Section  6.0)  and 
possibly  to  the  difference  in  the  pulse  train  envelopes.  Both  sets  are  in 
qualitative  agreement  with  theory  (21)  for  pulses  Gaussian  in  time  and 
space.  Although  the  theory  was  for  a  single  pulse,  the  pulse  train  can  be 
viewed  as  a  sum  of  separate  contributions  which  would  closely  approximate 
a  single  Gaussian  if  merged  together.  The  angular  range  of  incident  light 
in  these  experiments  was  less  than  0.1°  (FWHM). 

The  narrowing  of  the  phase-matching  curve  at  high  conversion  is  a 
factor  to  consider  in  the  design  of  practical  systems  from  two  points  of  view. 
Firstly,  the  effects  of  a  temperature  rise  following  crystal  absorption  become 
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Fig.  13  (a)  Low  conversion  phase  matching  curve,  (b)  68%  peak  energy 
conversion  narrowed  phase  matching  curve,  taken  with  3  psec  pulse 
trains. 
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Fig.  14  56%  peak  energy  conversion  narrowed  phase  matching  curve, 

contrasted  wath  low  conversion  case.  Taken  with  0.3  psec  pulse  trains. 


more  critical  and  secondly,  the  angular  acceptance  for  efficient  doubling  is 
restricted.  In  the  present  case,  with  a  3.5  cm  crystal  the  achievement  of 
>  50%  conversion  requires  that  external  full  angvilar  width  is  less  than 
about  0.13°.  Crystal  alignment  also  has  to  be  more  accurate. 
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6.0  HABMONIC  GENERATION  WITH  SHORT  PULSE  TRAINS 


Short  pulse  trains  obtained  with  the  grating  tuned  oscillator 
(Section  3.2)  were  employed  for  a  detailed  comparison  of  second  harmonic 
generation  with  theory.  The  duration  of  the  short  pulses  was  determined  to 
be  0.90  ±  0.05  nsec  by  autocorrelation  (Section  4.3  (b)).  The  pulse  trains  had 
a  full  width  at  half  maximum  of  0.30  psec  and  contained  about  10  pulses 
spaced  by  40  nsec.  The  maximum  energy  that  could  be  obtained  after 
amplifying  one  of  these  short  pulse  trains  was  143  mJ. 

Two  configurations  were  employed,  as  shown  in  Figure  15.  In  each 
case  the  transverse  profile  of  the  pump  beam  was  determined  at  the  crystal 
location  by  measuring  the  fraction  of  the  total  energy  transmitted  through  a 
range  of  different  iris  diameters.  This  data  was  least  squares  fitted  to  the 
transmission  function  anticipated  for  a  Gaussian  beam 

-2a^ 

T(a)=(l-e  w'  ) 

where  a  =  radius  of  iris  and  w  =  (l/e^  intensity)  radius  of  beam  (Section  2.3). 
Figure  16  shows  four  measured  profiles.  Parts  (a)  and  (b)  of  that  Figure 
refer  to  the  configurations  in  Figure  15  (a)  and  15  (b)  respectively,  although 
the  profile  for  Figure  16  (b)  was  taken  with  the  amplifier  turned  off.  The 
profiles  are  very  close  to  Gaussian,  with  l/e^  intensity  spot  parameters 
w  =  0.206  cm  and  w  =  0.220  cm  in  Figures  16  (a)  and  (b),  respectively.  When 
the  amplifier  is  turned  on,  the  near  field  profile  (configuration  of 
Figure  15  (b))  becomes  distorted,  as  shown  in  Figure  16  (c)  (solid  circles). 

For  comparison,  Figure  16  (b)  is  shown  again  in  part  (c)  (open  circles). 

This  distortion  is  due  to  the  amplifier  saturation  being  less  for  off-axis  parts 
of  the  beam,  leading  to  their  relative  increase.  Although  this  configuration 
gave  the  highest  fundamental  energy  (143  mJ)  the  presence  of  a  relatively 
large  low  intensity  "halo"  gave  lower  harmonic  generation  efficiency  than 
was  obtained  with  the  better  profiles  of  the  other  configuration 
(Figure  15  (a)). 

For  the  Figure  15  (a)  configuration,  examples  of  the  fundamental  and 
second  harmonic  pulse  trains  at  70%  energy  conversion  are  shown  in 
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Fig.  15  (a)  Focussed  doubling  configuration  with  5  m  lens;  (b)  Doubling 
configuration  with  direct  convergence  through  amplifier. 
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Figure  17.  The  depleted  fundamental  is  also  shown  in  Figure  17  (c),  not 
exactly  to  scale  because  of  a  changed  set-up. 

In  all  cases  the  second  harmonic  was  separated  from  the 
fundamental  by  passage  at  near-normal  incidence  through  an  imcoated 
1  cm  thickness  of  MgF2.  It  was  verified  using  a  spectrophotometer  that  this 
transmitted  less  than  0.5%  at  9.55  pm  and  about  92%  at  4.775  pm.  A  MgF2 
correction  factor  of  1.08  was  applied  to  all  the  second  harmonic  energies, 
which  were  measured  on  the  same  calorimeter  as  the  fundamental 
energies.  Energy  data  for  the  Figure  15  (a)  configuration  are  shown  in 
Figure  18. 

The  second  harmonic  conversion  formulae  of  Section  2.3  were 
applied  to  the  present  short  pulse  trains  assuming  that  the  pulse  train  itself 
had  a  Gaussian  envelope  in  time.  All  the  accurately  known  beam 
parameters  were  employed,  surface  entry  and  exit  losses  were  accounted 
for,  and  the  only  variable  was  the  nonlinear  coefficient,  dae,  of  AgGaSe2. 

We  show  three  values  of  dae  in  Figure  18,  to  illustrate  the  probable  range  of 
values  that  fit  the  present  data.  The  solid  line  and  two  dashed  lines  refer  to 
dae  =  31  pm  V-l  and  31  +  4  pm  V-l,  respectively.  Other  literature  values  are 
listed  in  Section  2.4  and  our  present  value  agrees  within  errors  with  two  out 
of  four  of  them,  tending  to  reinforce  the  lower  literature  values. 

On  the  beam  axis  the  peak  intensity  at  the  central  point  in  a  140  mJ 
pulse  train  is  160  MW  cm‘2.  In  the  above  configurations,  the  peak  fluence 
is  2.0  J  cm'2.  More  detail  is  given  on  the  surface  optical  damage  threshold 
in  Section  6.2.  At  the  fluence  of  2.0  J  cm'2,  there  was  visible  surface 
damage  within  three  shots.  However,  at  slightly  lesser  fluences  hundreds 
of  shots  could  be  taken  without  damage. 

6.1  Harmonic  Generation  with  Long  Pulse  Trains 

In  order  to  know  whether  the  surface  damage  threshold  is  related 
more  to  intensity,  or  more  to  fluence,  another  doubling  experiment  was  run 
with  a  much  longer  pulse  train  at  lesser  peak  intensity.  The  3  psec  pulse 
trains  shown  in  Figure  6  were  employed  in  the  optical  configuration  of 
Figure  15  (b).  The  beam  profile  at  the  crystal  is  shown  in  Figure  16  (d). 

This  profile  is  very  close  to  Gaussian,  showing  a  l/e^  (intensity)  radius 
w  =  0.188  cm.  The  difference  from  Figure  16  (c)  is  apparently  due  to  the  use 
of  a  larger  mode-selection  aperture  for  the  generation  of  long  pulse  trains 
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Fig.  17  Fundamental  and  harmonic  pulse  trains  at  70%  energy 
conversion,  (a)  9.55  pm  intensity;  (b)  4.775  pm  intensity;  (c)  depleted 
9.55  pm  beam  (not  exactly  to  scale). 
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Fig.  18  Calculated  and  experimental  second  harmonic  energy.  Solid 
curv’e  d36  =  31  pm  V'T  Dashed  curves  das  =  31  ±  4  pm  V’L 
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in  the  oscillator  (Section  3.2).  The  oscillator  spot  size  is  larger  and  hence 
the  far-field  spot  is  smaller. 

The  following  table  summarizes  the  harmonic  conversion  efficiency 
for  long  pulse  trains  at  varying  fluence. 


Fundamental 
Elnergy  (mJ) 

Harmonic 

Elnergy 

(mJ) 

Efficiency  (%) 

Fluence  (number  of  pulses) 

J  cm'2 

103.4  ±  2 

50.2  ±  3 

48.5  ±4 

1.86(10) 

156  ±  12 

114.6  ±5 

73.5  ±7 

2.80  (6)  occasional 
surface  spark 

200.7  ±  5 

157 

78 

3.6(1)  1  mm  diameter  of 
coating  damaged 

140 

(generation  of 
Phase  matching 

curve  at  62-68% 

conversion) 

2.5  (80) 

The  highest  harmonic  conversion  efficiency  of  78%  was  observed  on 
only  one  pulse  because  there  was  simultaneous  surface  damage  at  the 
entry  face  of  the  crystal.  Conversion  efficiency  of  62-68%  was  obtained 
during  a  phase-matching  measurement  at  a  fluence  of  2.5  J  cm*^  which  the 
surface  survived  for  at  least  80  pulses.  At  the  intermediate  fluence  of 
2.8  J  cm'2  there  were  occasional  small  visible  surface  sparks  on  the  crystal, 
which  however  did  not  lead  to  permanent  damage  because  the  same 
location  was  used  for  the  (subsequent)  phase  matching  measurement 
without  further  sparks. 

The  conversion  efficiencies  recorded  in  this  long  train  experiment 
were  just  as  high  as  those  for  the  short  train  experiment,  in  spite  of  the  use 
of  much  lower  peak  pulse  intensities.  One  reason  for  this  is  that  because  of 
the  profile  of  the  pulse  train  a  larger  fraction  of  the  pulses  are  of  sufficiently 
high  intensity  to  give  efficient  conversion.  Using  a  result  from  Section  2.3 
the  peak  Intensity  lo  in  an  individual  pulse  can  be  expressed  in  terms  of  the 
pulse  energy  Ep  by 
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,  _  2Ep  Iin2 

K 

where  A  is  the  (HWHM)  measure  of  pulse  duration  and  w  is  the  l/e^ 
intensity  radius  of  the  beam.  For  the  present  3  nsec  pulse  train  (75  pulses 
spaced  by  40  nsec)  the  average  peak  intensity  at  140  mJ  total  enei^, 
assuming  A  =  0.45  nsec  as  before  and  with  w  =  0.188  cm,  is 

lo  =  33MWcm"2 

From  Section  2.2  this  corresponds  to  a  peak  drive  parameter 

Vn^=CL7i^  =  2.04  using  d36=pmV“* 

and  an  average  conversion  efficiency  (using  Figure  1)  of  58%.  This  is  to  be 
compared  with  an  observed  efficiency  of  62-68%  at  the  same  fundamental 
energy.  The  apparent  discrepancy  could  be  due  to  a  reduction  in  pulse 
duration  through  the  pulse  train  because  of  continued  modulation 
(discussed  in  Section  4.2).  Although  an  autocorrelation  experiment  was  not 
performed  with  the  longer  pulse  trains,  there  is  evidence  from  the  direct 
detector/oscilloscope  measurement  that  the  pulse  duration  does  decrease 
somewhat  during  the  cause  of  the  pulse  train.  A  reduction  down  to  0.7  nsec 
(FWHM)  would  have  been  sufficient  to  give  the  observed  harmonic 
generation  efficiency,  but  this  lies  below  the  detector/oscilloscope  resolution. 

6J2  Damage  Hireshold  and  Fluence  Limitation 

Although  the  damage  threshold  for  AgGaSe2  is  often  quoted  as  an 
intensity  threshold  (typically  (22,14,23,15)  jn  the  range  10-30  MW  cm’2)  the 
present  experiments  with  short  pulse  trains  at  up  to  160  MW  cm'2  tend  to 
indicate  that  fluence  is  more  important  than  intensity  for  mode-locked 
pulse  trains  of  less  than  1  nsec  pulse  duration. 

The  harmonic  generation  experiments  were  performed  at  steadily 
increasing  fluence  until  visible  surface  damage  occurred.  Both  short 
(0.35  psec)  and  long  (3  psec)  pulse  trains  were  used,  in  order  to  differentiate 
between  the  effects  of  intensity  and  fluence.  Typical  pulse  train  envelopes  of 
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each  type  are  illustrated  in  Section  3.2.  In  each  experiment  the  transverse 
beam  profile  was  characterized,  as  discussed  above  in  Section  6.0,  and  the 
pulse  duration  was  measured  to  be  1  ±  0.2  nsec  in  general,  throughout  the 
measurements.  The  damage  observations  are  summarized  in  Table  2.  For 
the  short  pulse  trains  the  threshold  for  damage  fluence  lay  between  1.45  J 
cm‘2  and  1.58  J  cm'2  whereas  for  the  longer  trains  it  lay  between  2.5  J  cm'2 
and  2.8  J  cm'2.  In  view  of  the  approximately  10  times  greater  peak  intensity 
in  the  short  pulse  train,  the  relatively  small  reduction  in  damage  fluence 
(0.6  times)  indicates  that  fluence,  rather  than  intensity,  is  the  governing 
parameter  in  the  surface  damage  of  AgGaSe2.  This  suggests  that  a  critical 
surface  temperature  has  to  be  reached  for  damage  to  occur,  presumably  the 
temperature  at  which  sufficient  conduction  electrons  are  generated  in  the 
coating  material  to  contribute  a  significant  additional  9.55  pm  absorption. 
Over  the  range  of  these  measurements,  the  damage  fluence  can  be 
empirically  represented  by 


Fq  =  Jem  ^ 

where  t  is  the  duration  of  the  pulse  train  in  psec. 

In  other  recent  work(24)  the  damage  threshold  of  coated  and  uncoated 
AgGaSe2  was  measured  at  2.1  pm  using  180  nsec  laser  pulses.  An  initial 
set  of  coated  samples  gave  thresholds  ranging  from  1.4  2.5  J  cm  2.  After  re¬ 
polishing,  the  surfaces  that  had  been  re-coated  (by  a  different  vendor)  gave 
thresholds  ranging  from  2.8  -  3.3  J  cm'2,  whereas  the  uncoated  surfaces 
ranged  from  1.4  -  2.1  J  cm’2.  The  above  data  is  not  directly  relevant  to  our 
own  experiment  in  that  the  wavelength  and  pulse  waveforms  were  very 
different.  However,  it  illustrates  that  the  precise  details  of  surface 
preparation  and  coating  can  be  important. 


Tahle2 


Duration  of 

Pulse  Train 

(psec) 

Elnergyin 
Train  (mJ) 

l/e2 

Intensity 
Radius  (cm) 

Fluence 
(J  cm’2) 

Intensity 

(MWcm-2) 

Damage 

Y/N 

(no.  of  shots) 

0.35 

150 

0.49 

0.4 

50 

0.35 

180 

0.40 

0.7 

110 

0.35 

70 

0.21 

1.0 

190 

0.35 

100 

0.21 

1.45 

270 

0.35 

140 

0.24 

1.58 

295 

Y(l) 

3 

350 

0.49 

0.95 

20 

3 

140 

0.19 

2.5 

30 

3 

100 

0.19 

2.8 

33 

Y(3) 

3 

200 

0.19 

3.6 

43 

Y(l) 

6^  Optical  Quality  of  Second  Harmonic  Beam 

Tne  use  of  a  perfect  single  crystal  of  AgGaSe2  (Section  5.0)  and  a 
diffraction  limited  TEMqq  fundamental  beam  should  have  ensured  a  perfect 
second  harmonic  beam  quality,  and  indeed  this  turned  out  to  be  the  case. 

Using  a  range  of  iris  diameters,  the  fractional  transmission  of  the 
fundamental  and  second  harmonic  beams  was  measured  at  three  different 
distances  beyond  the  crystal:  0.18  m,  5.6  m  and  18.66  m.  These 
corresponded  to  the  near,  intermediate  and  far  fields  for  the  9.55  pm  or 
4.775  pm  radiation.  Typical  data  is  shown  in  Figure  19.  For  these 
measurements  the  hybrid  pulse  train  of  Figure  6  (a)  was  employed,  having 
60%  energy  conversion  efficiency  in  the  leading  spike,  less  than  10% 
conversion  in  the  tail,  and  approximately  30%  overall  energy  conversion. 

The  measured  l/e^  intensity  spot  parameters  are  summarized  in 
Table  (3). 

The  near  field  fundamental  spot  parameter  was  calculated  from 
Gaussian  beam  resonator  and  propagation  theory.^^S)  The  near  field 
second  harmonic  spot  was  calculated  in  the  approximation  that  walk-off 
distances  were  small  compared  to  the  fundamental  beam  diameter,  that  is, 
by  considering  the  conversion  of  a  series  of  pulses  Gaussian  in  time  but  like 
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a  plane  wave  in  space,  with  amplitude  appropriate  to  the  radial  location. 
The  necessary  formalism  is  discussed  in  Section  2.3. 

The  intermediate  and  far  field  fimdamental  spot  parameters  were 
calculated  from  resonator  theory,  plus  propagation.  In  the  far  field  the 
measured  spot  parameter  agreed  with  theory  to  within  4%.  For  the  second 
harmonic,  the  intermediate  and  far  field  spot  parameters  were  calculated 
from  the  measured  near  field  spot  parameter,  together  with  the 
assumption  that  the  phase  front  of  the  second  harmonic  had  the  same  (very 
slight)  curvature  as  the  fundamental.  Again,  the  agreement  was  very 
good,  although  in  the  far  field  the  measured  spot  was  11%  larger  than 
expected  on  the  basis  of  Gaussian  propagation  theory,  possibly  due  to  the 
small  fraction  of  the  harmonic  (less  than  16%)  that  originated  from  the  tail 
portion  of  the  hybrid  pulse,  and  would  have  been  generated  with  a  slightly 
smaller  spot  size.  The  far  field  spot  contained  90%  ±  10%  of  the  near  field 
energy. 
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Table  3  Summary  of  Beam  Quality  Data 


Distance  From  Crystal  (m) 
(Frequency) 

Measured  Spot  Parameter 
(cm) 

Calculated  Spot  Parameter 

0.18  (v) 

0.49 

0.50 

From  Resonator  Theory 

0.18  (2v) 

0.37 

-0.37 

From  SHG  Theory 

5.6  (V) 

0.72 

0.65 

Propagation  Theory 

5.6  (2v) 

0.51 

0.48 

Propagation  Theory* 

18.66  (v) 

1.45 

1.39 

Propagation  Theory 

18.66  (2v) 

1.08 

0.95 

Propagation  Theory* 

*From  measured  0,18  m  2v  spot  size. 
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7.0  THERMAL  CONDUCTIVITY  OF  AgGaSe2 

There  are  residual  absorptions  in  AgGaSe2  associated  with  either 
small  inclusions,  infrared  overtones,  or  band  edge  sites,  which  cause 
crystal  heating  under  high  average  optical  power  conditions.  The  crystal 
thermal  conductivity,  which  determines  the  transverse  (or  longitudinal) 
crystal  scale  size  that  can  be  used,  had  not  apparently  been  measured  prior 
to  the  present  work. 

A  simple  experiment  was  devised  to  measure  the  (relatively  low) 
thermal  conductivity.  A  mass  of  copper  (99.8  g)  was  contacted  via  the 
AgGaSe2  crystal  (10  mm  x  10  mm  x  35  mm)  to  a  temperature-controlled 
stage  as  shown  in  Figure  20.  Thermal  contact  was  assured  by  using  thin 
films  of  heat  sink  compound  between  the  aluminum  stage,  the  crystal  and 
the  copper.  Light  foamed  plastic  was  used  all  around  for  thermal 
insulation  and  the  temperatures  of  the  stage  and  the  copper  block  were 
monitored  by  digital  thermocouples. 

The  measurement  method  involved  a  period  of  stabilization  followed 
by  a  change  of  the  stage  temperature  to  a  new  value,  either  above  or  below 
the  onginal  value  by  10  to  15°C.  The  temperature  of  the  copper  block  was 
then  followed  to  its  new  equilibrium,  achieved,  in  the  case  of  perfect 
insulation,  only  due  to  heat  flow  through  the  crystal.  The  temperature 
difforcnce  between  the  copper  block  and  the  stage  decays  according  to 


ST  -  ATge - 


area  x  thermal  conductivity  V 
heat  capacity  x  crystal  width  j 


where  ATo  is  the  initial  temperature  difference.  Here,  the  area  of  the 
crystal  is  3.5  cm^  and  the  heat  capacity  is  that  of  the  copper  (38.44  J/°C)  plus 
a  contribution  from  the  crystal  itself.  The  heat  capacity  of  AgGaSe2  is^^G) 

1.8  J  cm'^/°C,  i.e.,  6.3  J/°C  for  the  entire  crystal.  Initially  there  is  an 
approximately  linear  temperature  gradient  between  the  stage  and  the 
copper  block  so  that  3.15  J/°C  is  the  effective  heat  capacity  of  the  crystal, 
where  temperature  change  on  the  block  side  of  the  crystal  is  concerned. 

The  heat  in  the  crystal  has  an  average  of  only  2/3  of  the  crystal  to  traverse, 
so  Its  effect  is  equivalent  to  only  2  x  3.15/3  =  2.1  J/°C  if  lumped  with  the  heat 
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P2538 


Fig.  20  Arrangement  for  the  measurement  of  thermal  conductivity. 
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capacity  of  the  copper.  The  combined  effective  heat  capacity  of  the  crystal 
and  the  copper  is  therefore  38.44  +  2.1  =  40.54  J/°C. 

Typical  temperature  decays  for  each  crystal  orientation  are  shown  in 
Figure  21.  In  one  case  the  conductivity  is  perpendicular  to  the  C-axis,  and 
in  the  other  it  is  at  (90-50.11)  =  39.89°C  to  the  C-axis,  i.e.,  mostly  parallel  to 
it.  There  is  no  significant  difference  between  the  temperature  decay  times 
for  either  of  these  orientations.  The  data  from  6  runs  are  therefore  lumped 
together  to  yield  an  average  experimental  decay  time  of  934.3  ±  26.4  sec. 
From  this,  the  conductivity  at  25°C  is 

K=  =0.0124 ±  0.0004 

3.5x934.3 

This  is  a  relatively  low  conductivity,  comparable  to  that  for  pyrex  or 
fused  silica.  Fortunately,  the  temperature  dependence  of  the  phase 
mismatch  for  doubling  to  4.775  pm  is  almost  zero  (Section  5.0),  so  that 
temperature  excursions  of  more  than  100°C  can  readily  be  tolerated. 
Further  discussion  of  the  use  of  AgGaSe2  for  high  average  power  frequency 
multiplication  is  given  in  the  following  section. 
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Fig.  21  Two  examples  of  measured  copper  block  temperature  decays  and 
their  fitting  to  an  exponential  function. 
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8.0  mOH  AVERAGE  POWER  DOUBLING 


In  this  section,  we  will  briefly  address  some  of  the  issues  relating  to 
high  average  power  doubling  in  AgGaSe2-  The  only  limitation  is  crystal 
heating  due  to  absorption  at  either  9.55  pm  or  4.775  pm.  For  the  typical  case 
in  which  the  length  of  the  crystal  greatly  exceeds  its  transverse  dimension 
a  radial  temperature  profile  is  set  up  between  the  axis  of  propagation  and 
the  outside  surface  of  the  crystal.  The  refractive  index  increases  with 
temperature  C15,20)  go  ^  positive  lens  is  formed.  At  sufficiently  high  internal 
temperature,  differential  thermal  expansion  causes  stresses  which  will 
ultimately  crack  the  crystal.  An  estimate  of  the  maximum  temperature 
differential  for  AgGaSe2  is  100-200°C,  derived  by  Eimerl  et  al.^27) 

In  order  to  estimate  the  temperature  rise,  consider  a  cylindrical 
crystal  of  radius  a  and  thermal  conductivity  K  subject  to  a  uniform  volume 
heating  rate  Ao  (W  cm'^).  If  the  length  of  the  cylinder  is  much  greater  than 
its  diameter  and  the  outer  cylindrical  wall  is  kept  at  temperature  To  then 
the  steady  state  temperature  profile  as  a  function  of  radius  is^^S) 

T(,)=T„+^ai£za 

In  reality  the  heating  rate  will  not  be  uniform  throughout  the  crystal 
because  the  beam  intensity  varies  radially  and  the  fundamental  and  second 
harmonic  absorption  coefficients  will  be  different,  leading  to  a  variation 
with  length.  Although  the  present  experiments  used  a  Gaussian  radial 
intensity  profile,  it  would  be  better  from  the  standpoint  of  efficient 
utilization  of  crystal  material  to  use  a  "super-Gaussian "  profile,  which 
begins  to  approximate  a  constant  volume  heating  rate.  As  it  happens,  the 
9.55  |xm  and  4.775  pm  absorptions  in  the  present  AgGaSe2  sample  are  both 
about  1.5%  per  cm,  so  that  a  rough  approximation  to  the  crystal 
temperature  rise  can  be  obtained  by  assuming  that  the  heat  input  is  also 
constant  along  the  crystal  length. 

For  the  above  temperature  profile  there  is  a  refractive  index  profile 

n(r)=no+a(T(r)-To) 
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where  n(r)  is  the  index  at  radius  r,  no  is  the  index  at  radius  a  and 
temperature  Tq,  and  a  is  the  coefficient  of  refractive  index  with 
temperature.  This  is  only  an  approximate  estimate  because  in  a  cylindrical 
geometry  there  are  stresses  leading  to  compression  of  the  central  material 
and  rarefaction  of  the  outer  shell,  resulting  in  a  modified  index  profile.  The 
expansion  of  AgGaSe2  is  anisotropic^29)  having  in  fact  a  negative  coefficient 
parallel  to  the  c-axis  and  a  positive  coefficient  perpendicular  to  the  c-axis. 

Considering  only  the  thermal  dependence  of  the  index  the  focal 
length  of  the  equivalent  lens  can  be  calculated.  For  a  thin  lens  of  radius  a 
with  one  surface  plane  and  the  other  of  radius  R  the  thickness  is  closely 
approximated  by 


If  it  is  made  of  a  material  of  refractive  index  p,  the  optical  thickness  at 
radius  r  is  t(r)  (p-1),  which  is  to  be  equated  to  Lc  [n(r)-noJ  from  the  thermal 
lens  in  a  crystal  of  length  Lc.  Using  the  above  expressions  and  the 
relationship  f  =  R/(p-l)  for  the  focal  length  f  of  the  thin  lens,  we  obtain 


for  the  focal  length  of  the  thermally  induced  lens  in  the  crystal.  In  terms  of 
the  absorption  coefficient  of  the  crystal,  Uq,  the  heat  input  is  Aq  =  loCx©  where 
lo  is  the  (assumed  uniform)  incoming  laser  intensity.  Then 


For  example,  with 
K  =  0.0124  Wcm-lK-I 

Lc  =  3.5  cm 

a  =  30  x  10-6  K-l 

lo  =  200  W  cm-2 

tto  =  0.015  cm- 1 

we  have  fc  =  79  cm. 


L=- 


2K 


Lpalna 


0^*0 
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r 


The  temperature  rise  at  the  center  of  the  crystal  relative  to  the 
(cylindrical)  edge  is 


T(o )  -  To  = 

°  4K  4K 

This  now  depends  on  the  crystal  radius  a-  For  the  above  example, 
with  a  =  0.5  cm,  the  central  temperature  rise  is  15°C. 

The  effect  of  temperature  upon  phase-matching  is  discussed  in 
Section  5.0  above.  As  shown  there,  the  doubling  process  is  completely 
insensitive  to  temperature,  at  least  over  the  first  50°C  above  room 
temperature.  It  appears  that  the  most  significant  phenomenon  at  high 
average  power  will  be  the  formation  of  a  positive  lens  in  the  crystal.  To  the 
extent  that  this  lensing  action  begins  to  induce  angular  deviations  in  the 
beam,  the  phase  matching  will  be  affected,  particularly  at  high  conversion 
efficiencies  (Section  5.1).  For  phase  matching  to  be  preserved,  the  angle  a/fc 
should  remain  less  than  about  0.2  degrees  (3.5  x  10'3  radius). 
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